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Abstract 

In the framework of the QCD string approach it is shown that the spin-averaged masses M{nL) 
of all low-lying light mesons are well described using the string tension a as the only parameter. 
The Regge slope a! L and the intercept ol(0) of the Regge L-trajectory for M(nL) are calculated 
analytically and turn out to be a' L = 0.80 GeV~ 2 (for L < 4) and ol(0) = —0.34, in good 
agreement with the experimental data: a' L = 0.81 ± 0.01 GeV~ 2 , aL ex p(0) = —0.30 ± 0.02. To 
obtain this strong agreement with the data the nonperturbative quark self-energy contributions to 
the meson masses must be taken into account, which appeared to be large and negative for small 
values of L, and are even for larger values of L important for a close fit. From the present analysis 
of the meson spectra the restriction a s < 0.40 on the strong coupling is required. 

PACS: 12.38.Lg, 12.39.Ki, 14.40. Cs, ll.10.Ef, 11. 10. St, 11.15.Tk 



I. INTRODUCTION 



The spectra of hadrons form an extremely important test ground for nonperturbative 
QCD. The scaling property of QCD tells us that in the end all characteristics of hadrons 
must depend on a single parameter, say o or Aqcd- Till now, all calculations of hadronic 
spectra with an accuracy comparable to the uncertainties in the experimental data, have 
relied on models that contain several, in some cases even many, parameters. Here we adopt 
the formalism that uses the QCD-string Hamiltonian which relies on only one parameter- 
the string tension a, while for light mesons the interactions derived from perturbative QCD, 
both static and spin-dependent, can be considered as a perturbation. The string tension can 
be extracted from experiment, in particular from the slope of the leading Regge trajectory, 
and in the present paper we use it to describe the orbital excitations of light mesons. 

The QCD string approach developed in recent years Q - @ starts from first principles, 
i.e. from the Euclidian QCD Lagrangian. In Ref. the relativistic Hamiltonian Hr for the 
light mesons with spinless quarks was derived under several verifiable assumptions. 

First, string excitations (hybrids) are not taken into account, since there is a large gap, 
~ 1 GeV, between a mesonic ground state and its gluonic excitation [[|. Without this 
approximation one obtains a matrix multichannel Hamiltonian ||. Therefore the ground 
states of the mesons with not too large orbital momentum [L < 5) can be treated in the 
closed channel approximation. 

Secondly, the relativistic Hamiltonian Hr used here is derived in the quenched approxima- 
tion, where creation of qq pairs (sea quarks) is neglected. The accuracy of this approximation 
for low-lying states is expected to be approximately 10% while high radial and orbital 
excitations can be strongly affected by qq pair creation|J. 

Thirdly, to derive the Hamiltonian only forward- in-time quark trajectories have been 
taken into account while backward trajectories of a quark (antiquark) were neglected. The 
accuracy of this approximation has been checked comparing meson spectra of the center-of- 
mass Hamiltonian and light-cone Hamiltonian H^c, since in the latter case backward- 
in-time trajectories do not contribute M. This comparison has shown that the differences in 
meson masses for these Hamiltonians is not larger than 10% for all mesons with exception 
of the pion, which receives a large contribution from backward-in-time trajectories corre- 
sponding to negative energy components and for them the formalism should be modified. 
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In the Hamiltonian used chiral effects are not taken into account. 

Under these assumptions and with the use of the Fock-Feynman-Schwinger (FFS) rep- 
resentation the Green's function G of a meson with a spinless quark and antiquark can be 
written in a gauge invariant way as a functional integral with the action A M: 

G = ds ds DzDz exp(— A), 
Jo Jo 

A = K + K + aS min , (1.1) 



where in the action (|1 . 1| ) the only approximation made is that the vacuum average over the 
Wilson loop (W(C)) is taken in the form of the minimal area law, viz 

(W(C)) = const exp(-cxS min ). (1.2) 

The accuracy of this approximation is determined by the condition 

i?>T 9 , (1.3) 

where R is the size (e.g. the r.m.s. radius) of the meson and T g is the gluonic (vacuum) 
correlation length which determines how the vacuum correlators decrease as a function of 
the separation r between the quark and the antiquark. The value of T g was calculated on 
the lattice || and in the quenched case was found equal to ~ 0.15 — 0.20 fm, i.e. it is 
much smaller than the r.m.s. radii R(nL) for light mesons having R(nL) > 0.8 fm. Note 
that the condition (|1.3|) is also valid for most excited heavy-light mesons and even in heavy 
quarkonia, where e.g. R(lP,cc) ~ R(2P,bb) ~ 0.65 fm. 

In the general case the approximate area law can be replaced by the exact qq interac- 
tion which contains a linear confining part for distances satisfying Eq. ( |1.3| ) but should be 
modified at smaller distances, being fully determined by a bilocal field correlator. In the 
action A (|1 . 1|) K and K are the quark (antiquark) kinetic energy term written in the FFS 



representation: 



K 



m 2 + - A z{r) 2 a 



dr, K= m 2 + -z{r)l dr, (1.4) 
Jo L 4 

where m is the current mass of a quark, r is the proper time introduced by Schwinger[^], 
z(t) and z(t) are the paths of the quark and antiquark. 

To define the Hamiltonian one can use the connection between the meson Green's function 
and the Hamiltonian H: 

OG/dT = -H R G, (1.5) 



where the Hamiltonian can be determined on any hypersurface, i.e. H can be derived in 
different frames. Here we shall use the Hamiltonian obtained in the cm. frame, while in 
Ref. ]7| the Hamiltonian was derived in the light-cone frame. 

In order to use the relation ( |1.5| ) in Euclidean space-time it is of great importance to go 
over from the proper time r to the actual time t = 04, of a quark (antiquark). Doing so the 
new quantity ji{t) is introduced: 

2/i(t) = dt/dr. (1.6) 

After perfoming the canonical quantization the variable /1 (being a canonical coordinate) 
will define the constituent mass of a quark. The last term in the action ( |1.1[ ) has the form 
of the Nambu-Goto string: 

S min = J* dt £ dp y / det^, (1.7) 

where g a b = d a w a dbW a with a,b = t,j3 and w a (t,j3) are the coordinates of the string world 
surface. In Refs. the string was approximated by a straight line connecting the path 
coordinates z a (t) and z a (t) and in this case 

w a (t, (3) = z a (t){3 + z a (t)(l -(3), < < 1. (1.8) 

Due to the presence of the square root in S m { n in Eq. ( |1.7| ) this term cannot be quantized 
and to get rid of it one can introduce two auxiliary fields, u(t, P) and 7](t, P) in the standard 
way as done in string theory [10 . 



By definition the introduction of the auxiliary fields is accompanied by the additional 
integration over Dfi, Du, Dt] in the functional integral defining the meson Green's function, 

G = J DfiDuDi]DR a Dr a exp(-A), (1.9) 

where "the center-of-mass" coordinate R a and "the relative" coordinate r a are introduced 
instead of the path coordinates z a and z a . 

As shown in Refs. the integrations over Dr\ and DR a can analytically be performed 
in the integral (|1.9|) and after that the Green's function has a simpler form, viz 

G = J D^DvDr exp{-A R ). (1.10) 

As the next step instead of performing the integration over Dfi and Dv in the integral 
( PTH ) one can use an equivalent procedure-to go over to the canonical quantization of the 
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Hamiltonian Hr which corresponds to the action Ar. This Hamiltonian in Minkowski space- 



time is obtained from in a standard way|2|, |10[ and is taken as a starting point in our 
analysis (see Eq. (|3.1|)). 

The quantization of has been performed in two cases: the quasiclassical quantization 
of Hn for large orbital momenta L, L 3> 1, in Ref. [[□]] and for not too large L in Refs. [0. 



Here we are interested only in orbital excitations with L < 4 and in this case the Hamiltonian 
can be presented as the sum 

H R = + AH stT , (1.11) 

where the term AH str is relatively small for L < 4 and can be considered as a perturbation. 
Then the problem reduces to the quantization of more simple Hamiltonian which can be 
easily done (see Sections |HI|and |V|). At the final stage we obtain the surprising result that 
after quantization the unperturbed Hamiltonian coincides with the Hamiltonian used 



in the relativistic potential model (RPM) [12], [13]. In this way (for the states with L < 4) 
the connection between the QCD string Hamiltonian Hjp and the RPM is established and 
we can also calculate the corrections which are absent in the RPM. 

The first correction A str , called the string correction, comes from the term AH str in 
Eq. ( |1.11| ); it is negative and varying from a value of about -50 MeV for L — 1 to about 
-150MeV) for L = 4. 

A second correction to the meson mass, Ase, is due to the spin (colour magnetic moment) 
interaction of a quark (antiquark) with the external (vacuum) field when the operator cr Q( gF Q( g 
is inserted in the Wilson loop. The form of this operator is 

° a p = \ (7«7/3 - ipla)- (1-12) 



This nonpertubative self-energy correction was analytically calculated in Ref. |T4| with the 
the use of the FFS representation for the quark Green's function. A SE is negative and has 
a rather large magnitude (of the order of -400 MeV to -300 MeV) for all states with L < 4, 
slightly decreasing with growing L. Owing to this correction the correct value of the Regge 
intercept was obtained. 

We use here the current quark mass m = and due to the procedure of canonical 
quantization the constituent mass of a quark can be defined in a rigorous way as the canonical 
coordinate [i which is equal to the quark kinetic energy. It is of interest to note that the final 



expression of the Hamiltonian does not contain the constituent mass at all. However, the 
constituent mass /i must be defined since it explicitly enters those terms in the Hamiltonian, 
like the spin-dependent and self-energy terms, which are considered as a perturbation. 

Here we consider in detail the spin-averaged meson masses, M(nL), or the centers of 
gravity of the nL-multiplets (i.e. neglect the hyperfine and fine-structure splittings) for which 
the physical picture is simpler and at the same time more universal since the parameters do 
not depend on spin and isospin. 

We concentrate mostly on the orbital excitations with n = for which experimental 
data exist for all ground states with L < 5. Then for the linear confining potential or 
all meson masses M(nL) can be expressed through a single parameter-the string tension 
a . The values of the slope and the intercept of the Regge L-trajectory (L < 4) will be 
calculated analytically: their numerical values are a' L = 0.80 GeV~ 2 (a = 0.18 GeV 2 ) and 
Q! l(0) = —0.34 turn out to be in very good agreement with the experimental numbers. 
From the Regge slope a restriction on the admissable values of the string tension follows: 
a = 0.18±0.005 GeV 2 for the pure linear potential and a = 0.19±0.01 GeV 2 if the Coulomb 
interaction is taken into account. 

The Coulomb contribution is mainly important for the IS" and the IP states having 
values in the range —200 to —100 MeV, and is considered here in a twofold way: from exact 
calculations with the linear plus Coulomb potential and also when the Coulomb interaction 
is considered as a perturbation; both considerations give very close results. For a = 0.19 
GeV 2 the QCD coupling is a s = 0.39 which is typical for heavy quarkonia, and from our 
analysis of the meson spectra the following restriction on the strong coupling a s < 0.42, 
is obtained. This number is in a good agreement with the two-loop value of the freezing 
coupling constant obtained in background field theory [15| . 



II. EXPERIMENTAL DATA 

The experimental numbers for the spin-averaged masses M(L), or the centers of gravity 
of the 1 L-multiplet (the ground states with n = 0), are presented in Table | and need some 
remarks. 

First, all members of the 1 3 Pj multiplet are supposed to be known: 02(1318), ai(1235), 
and the ao(980) too, are considered to form the 1 3 Pq multiplet with M(1P) = 1252 MeV. 
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a 6 (2.45) qd /U2.51) 




FIG. 1: The Regge L-trajectory for the light mesons. 



Similarly, for the fj(lP) mesons the spin-averaged mass is 1245 MeV fl~6| . 

Second, in the case of L = 2, the fine-structure splittings of the 1 D-w&ve mesons are 
supposed to be suppressed as compared to the P-wave states |T7J. As a result, for all 
members of the 1 Dj-multiplet, e.g. the ^3(1. 69) and 7^(1.67), their masses are very close to 
each other and one can expect that the true value of M(1D) lies between these two values. 
The same would be valid for the isoscalar mesons, if they were not mixed with other hadronic 
states, and just this situation is observed in experiment where the masses of the a>3(1.67) 



and the cj(1.65) have values close to the corresponding isovector mesons ||T6 



Due to the suppression of the matrix elements (m.e.) like < 1/r 3 > the spin splittings 
for the higher orbital excitations like IF, 1G, etc. should be even smaller than for the ID 
mesons. Therefore the masses of the 04(2.01) and the /4(2.03) are supposed to be close to 
M(1F) as well as the masses of the ps(2.30) with L = 4, and ag(2.45) and /g(2.47) with 
L = 5, lie close to their centers of gravity. The masses of all orbital excitations (n = 0) can 
be nicely described by the Regge L-trajectory (see Fig.l): 



M 2 (L) = (1.23 ± 0.02) L + 0.37 ± 0.02 (GeV 2 



(2.1) 
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or 

L = 0.81M 2 (L) — 0.30 (2.2) 
with the following Regge slope and intercept: 

a' Lexp = 0.81 ± 0.01 (GeV 2 ) and a Lexp (0) = -0.30 ± 0.02 (L < 4). (2.3) 

and both values have a small experimental error. 
Note that for the leading p-trajectory: 

j = a'jM 2 {J) +0.48 (2.4) 

the slope a'j = 0.88 GeV -2 and the intercept «j(0) = 0.48 are larger since their values 
depend on the spin contributions. On the contrary, the L-trajectory is a universal one and 
in the approximation of closed channels it is the same for isovector and isosinglet mesons. 

In our paper the meson masses, Regge slope and Regge intercept will be calculated 
analytically in the framework of the QCD string approach. 



III. RELATIVISTIC HAMILTON! AN 

We start with the Green's function Eq. ( |1.10j ) which is obtained after performing two 
integrations DrjDR in the functional integral, so in the Green's function still three integra- 
tions left: over the two auxiliary fields fi(t) and v{t) and also over the separation r between a 
quark and an antiquark, r a = z a (t) —z a (t). The explicit expression of this action was ob- 
tained in Refs. and the dependence of Ar on \x and v turns out to be rather complicated. 
Therefore it is more simple and convenient to go over to the equivalent procedure-canonical 
quantization of the Hamiltonian Hr which corresponds to the action in Eq. (|1.10|) : 



p 2 , + m 2 

h r = — r~^~ + ^ r 

/IT 



+ — 



Mr) + 2 I' dPv{P){P-\Y 
Jo 2 



By definition the field operator fi{t) is 
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where t is the actual time. 

In Eq. m is the current quark mass, which for a light quark (antiquark) will be 

taken equal to zero; L is the angular orbital momentum, L = f x p, and the operator 
Pr = (p ' r) 2 /(r 2 ). Canonical quantization of H R has not been done in general, but 
was performed in two cases: first when the orbital momentum L ^> 1 (the quasiclassical 
quantization in Ref. flTT| ) and for relatively small L < 4 which will be considered here. The 
constant o determining the nonperturbative potential is the string tension. 

It is convenient to rewrite H R as a sum of two terms: 

H R = H$ + A# str , (3.3) 
with the "unperturbed" Hamiltonian defined by 

where in we have included the term L 2 / (fir 2 ) and subtracted the same term to give the 
string correction AH stT 

L 1 



AH str = ~ 



r 2 [/x(t) ^ + 2tid{3v{{3){(3-\) 2 



2 $ df3v{(3){P ~ \f (35) 



fir 2 fi + 2f 1 df3v(f3)([3 



1)2- 
2> 



If L is not large, then the term AH stI appears to be relatively small and can be considered 
as a correction to the Hamiltonian H^' [11] but for large L the representation of as the 



sum Eq. ( |3.3j ) is of no use, since in this case both terms are equally important. Note that 
to get the expression (|3.4| ) one needs the following definition 

L 2 



P 2 =Pr + ^2- (3 - 0! 



The simplest Hamiltonian H with L = is a special case of H R (or H# ) with p 2 replaced 
by p 2 . 

Since AH str is considered as a perturbation, the canonical quantization needs to be per- 
formed only with the Hamiltonian which can be easily done, since the latter can be 
presented as the sum of the kinetic energy operator ifkin which depends on \x and V^ t which 
depends on v. 

^ 1) = F kin (/i) + t> pot (z/). (3.7) 



The variables /x and v are the canonical coordinates while the momenta ir^ and n v canonically 
conjugated to /i and v turn out to be equal to zero, since (as well as H-r) does not 
depend on the derivatives fi and v. Hence it follows that 

^ = dH^/d/j = 0, tt u = dH^/dv = 0. (3.8) 

Thus instead of calculating the Green's function as the functional integral Eq. ( |1.1(J| ) and 
subsequently deriving the Hamiltonian H# , one can obtain this Hamiltonian from 
Eq. fl3.4| ) with the use of the extremum conditions Eq. (|3.8j ). The equivalence of these 
two procedures was explicitly illustrated in the third paper of Ref. [@ for the kinetic part in 
Eq. 0: 

HUV) = P \ m + (3-9) 
For this Hamiltonian the explicitly calculated Green's function for the free particle is 

G = J Dz(t)Dpexp[i j\p '■ 'z - y 'p 2 + m 2 )] (3.10) 

which is just the canonical representation of the Green's function with the free Hamiltonian 



H = ^fp 2 + m 2 (3.11) 

IV. THE EXTREMAL VALUES OF THE OPERATORS /j AND v 

To understand the physical meaning of the auxiliary fields /i(t) and v(i) let us find their 
extremal values. First, in the Hamiltonian we determine the variable v{f3) from the 
extremum conditions ( |3.8| : 



- , 0, — -f- = 0. (4.1) 

Then one finds that v{(3), which is an operator in general, does not depend on the string 
parameter (3 and is equal to 

MP) = °r, (4-2) 

i.e. is actually the energy density along the string. 

With the use of Eq. (|4.2|) the Hamiltonian H^p reduces to a simpler operator: 

H® = + p{t) + ar, (4.3) 
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where fi(t) is still an operator in the Hamiltonian formalism. Its extremum can be found 
from the second extremum condition (TO) 



/i(t) = v /p 2 + m 2 ,for^ 1) , f j lQ (t)=y/p^+m i ,ioTH , (4.4) 

i.e. the extremal value of /i is one half the kinetic energy operator. Note that after canonical 
quantization the 'coordinate' fi is already independent of time. Substituting it into the 
Hamiltonian one obtains 



= + m 2 + or (4.5) 

giving rise to an eigenvalue equation that is identical to the spinless Salpeter equation (SSE) 
with a linear potential 

H^ip(nL) = M (nl)ip(nL). (4.6) 



This equation has been used in the RPM for many years fl"2|, |i~3| and the only difference is 
that in the Hamiltonian in Eq. ( |4.6[ ) we must use the current quark mass m since just 
the current mass enters the meson Green's function (m = for the light mesons) while in 
the RPM m ^ is usually used, e.g. in Ref. @ m = (m u + m d ) /2 = 220 MeV. 



V. THE CONSTITUENT QUARK MASS 

Although the constituent mass /i is not explicitly present in it enters many important 
physical characteristics like the spin splittings and magnetic moments, and also in the string 
and self-energy corrections, therefore it must not be left in as an operator. The simplest 
way to solve this is to define /j,q as the expectation value of one half the quark kinetic energy 
operator Eq. (|4.4j) , i.e., 

HainL) = (yp2 + m 2 ) nL . (5.1) 

Note that the eigenvalues Mo(nL) in Eq. Q4.5Q for the linear potential or are connected with 
Ho as follows 

M (nL) = 4^o (nL). (5.2) 

The values of //o can be expressed through a single parameter-the string tension o and the 
universal numbers a(nL) given by 

fi (nL) = y/oa(nL). (5.3) 
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This relation is a manifestation of the scaling property of the SSE in the case m = 0. 

Another definition of the constituent mass, denoted by fto, was used in Refs. ]!], [3] in the 
so called "einbein approximation" (EA) where the second extremum condition in Eq. ( |3.6| ) 
is written not for the operator but for the eigenvalues M (nL). A priori it is not clear 
whether in both definitions the extremal values ^{nL) and fto(nL) coincide or not, therefore 
let us compare them. In the EA Eq. ( [4.6[ ) is rewritten as 



p 2 + m 2 

+ ar 



i) = e{nL)ip, (5.4) 



with 

e(nL) = f — J A(nL) } (5.5) 

i.e., it reduces to the Airy equation with e(nL) = M (nL) — ft and the quantities A(nL) in 
Eq. ( |5.5|) are the zeros of the Airy function. The constituent mass jlo(nL) is now determined 
by the condition 

^M + i = , (m = 0). (5.6) 
aft 

Then from Eqs. flOj ^6|) one obtains that 

n \3/4 

Ho{nL) = ^/a ( -A(nL)\ = y/aa{nL). (5.7) 

To compare fio(nL) and fto(nL) one can use the numbers presented in Appendix [A| (see 
Tables [VTTI and |Vlll[ ) from which the corresponding universal numbers a(nL) and d(nL) can 



be determined. 

The largest difference between n Q (nL) and fto(nL) was found for S waves and is increasing 
with growing radial quantum number n from 5% for the 15* state to 7% for the 55* state. 
However, this difference is falling with increasing L, being only 1.7% for L — 5 (n = 0). So, 
fio and fto are numerically very close. In contrast to the eigenvalues Mo(nL) for the Salpeter 
and Airy equations a large difference is found between some matrix elements (m.e.) like 
< 1/r 3 > (for any L ^ 0) which define the fine-structure splittings. This difference can be 
as large as 30-50% in some cases (see Tables VII , |VIII| ). Moreover, while for the SSE these 



m.e. are growing, they are slightly decreasing for the Airy equation. It is worth to notice 
that these differences between the m.e. would be much larger if a fixed constituent mass, as 
in potential models, would be used. 
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The reason behind such discrepancies may be connected with the different asymptotic 
behavor of the wave functions(w.f.). For the SSE Eq. Q4.5Q it falls as exp(— \for) ||18| while 
for the Airy Eq. ( |5.4j ) the w.f. decreases as exp(— jl a r 3 ^ 2 ) . Therefore the definition (|4.6|) 
of the constituent quark mass as well as the calculations of the m.e. with the use of the 
unperturbed Hamiltonian H^/ has to be considered as preferable compared to the EA. 

Note a useful relation between the m.e.: 

(err) = 2fi (nL) (5.8) 

and 

(1/r) = MVp)nL, (5.9) 

where (1/p) is independent of a but does depend on the quantum numbers. 

It is worthwile to discuss some common features and differences between the QCD string 
approach used here and the RPM which was an essential step in our understanding of 
hadronic spectra. 

First of all we have shown that the Hamiltonian used in the RPM is more than a model 
one, as it can be deduced from the meson Green's function in QCD for not too large angular 
momenta L, assuming that the string correction coming from the part AH str Eq. (|3.5|) is 
neglected or considered clS db perturbation. 

However, in the QCD string approach the mass m must be the current quark mass, 
since just the current mass enters the meson Green's function in the FFS representation, in 
particular for light quarks m = and for the strange quark m s = 140 — 160 MeV are taken 



@, while in the RPM m = (m s + m d )/2 = 220 MeV, m s = 465 MeV are taken in p . 

Since in the QCD string approach the Hamiltonian as well as the string and self-energy 
corrections are calculated with the use of just the same FFS representation the whole picture 
is simplified and the spin-averaged mass can be expressed through the only parameter-the 
string tension a. 

The constituent mass of a quark does not enter the final form of the Hamiltonian Eq. ( |4.5|) 
and the notion of the constituent mass appears to be necessary only when one takes into 
account the spin-dependent interaction and also the string and self-energy corrections to the 
meson mass. In the FFS representation all these terms are inversely proportional to just 
the same auxiliary field fi |l[] which is strictly detremined from the extremum condition due 
to the procedure of canonical quantization and appears to be the kinetic-energy operator. 
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TABLE I: The masses Mq(L) and Mg(approx) = 8L + 3tt(t for the ground states (n = 0) with 
L < 6 (a = 0.18 GeV 2 ). 



L 





1 


2 


3 


4 


5 


6 


M 2 (nL) 


1.7940 


3.2126 


4.6431 


6.0777 


7.5142 


8.9518 


10.3900 


8La + 3ira 


1.696 


3.1365 


4.5765 


6.0165 


7.4565 


8.8965 


10.3365 


difference 


5.4% 


2.3% 


1.4% 


1.0% 


0.77% 


0.62% 


0.5% 



However,to calculate these corrections the variable jj must not be used as an operator and we 
have defined the constituent mass as the expection value of the quark kinetic energy. This 
definition of the constituent mass is in accord with another one-the variational definition of 
the constituent mass used in Ref. 

It is important that the constituent mass in our case depends on the quantum numbers 
and increases with growing L and n r . 

Finally, instead of the string and self-energy corrections considered in next Sections (they 
are negative) in the RPM a universal negative constant is introduced. 

VI. THE STRING CORRECTION AND THE SLOPE OF THE REGGE TRAJEC- 
TORY 

It is known that for the Salpeter equation ( |4.6| ) (or for the unperturbed Hamiltonian 
H^) the squared masses Mq{tiL) can be approximated (with an accuracy of about 1% for 
L + 0) by the "string formula" 0, 

M 2 (approx) = 8a L + 4vra(n + 3/4). (6.1) 

The exact values of Mq{hL) together with those of M^approx) (L < 6 , n = 0) are given 
in Table | from which one can see that the differences between them are indeed < 1% for 
L > 2. 

As is clear from the approximation (|6.1|) , the slope of the Regge trajectory for the SSE is 
(8<t) _1 , i.e. a' L = 0.69 GeV~ 2 for a = 0.18 GeV 2 , which is 17% smaller than the experimental 
number Eq. ( |2.3|) , a' L = 0.81 ± 0.01 GeV -2 . Note that the string corrections which come 
from the term Eq. Q3.5|) are also proportional to L and therefore affect the Regge slope. The 
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TABLE II: The string corrections A s t r in MeV and the mass Mq(L) in GeV, for the ground states 
(L<6). 



L 


1 


2 


3 


4 


5 


6 


M (L) 


1.7924 


2.1549 


2.4653 


2.7412 


2.9920 


3.2234 


A str (L) 


-52.9 


-86.9 


-113.0 


- 132.7 


-153.7 


- 170.7 


A str (asym) a 








-142.4 


-182.9 


-219.1 



a ) For L < 3 the asymptotic formula Eq. ( |6.12| ) is not applicable. 

situation appears to be different in two domains: L < 4 and L > 5 respectively, and we 
consider them separately. 

A. Case A. L < 4 



By the definition (|3.5| ) AH stI gives a negative correction to the eigenvalues M (nL); its 
magnitude turns out to be relatively small, ~ —100 MeV, and therefore this term can be 
considered perturbation |TT| . 

aL(L + 1) / 1 \ 



A stI (nL) = (AH stI ) 



(6.2) 



fio(nL) \r(6fi + ar) / ' 

In Eq. ( |6.2|) we have used that the integral Jq df3(f3 — 1/2) 2 is equal to 1/12 and the operators 
v and fi were replaced by their extremal values Eq. (|4.2|) and Eq. ( |4.4|) . The factor in brackets 
can also be approximated (with an accuracy better than 3%) replacing or by < or >. Then 
the string correction is 

A , a(nL) = jmm , (6 .3) 



M 6 /^) + (<"*)) ' 

Due to the relations ( |5.2| ) and ( |5.8| ) for the linear potential the correction A str becomes 

2a 3 / 2 (l/p)L(L + l) 



oL{L + l){l/r) 



*str 



Ml 



(6.4) 



Note that in Eq. ( |6.4|) the m.e. < 1/p > VX + I is almost constant, varying from 0.787 
for L = 1 to 0.741 for L = 4 (see Table [VI I|) . The values of A str (using the m.e. < 1/r > 
from Tables VIII and [IX| ) are given in Table |T|. 



For comparison in Table [TJ the string corrections valid for large L (L > 5) (see the 
asymptotic string correction formula Eq. ( |6.13 )) are also given. 

Now one can analytically calculate the Regge slope for the "corrected " mass: 
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TABLE III: The squared masses M 2 = (M + A str ) 2 in GeV 2 for the ground states (L < 6, a = 0.18 
GeV 2 ). 



L 


1 


2 


3 


4 


5 


6 


(M + M str ) 2 


3.026 


4.277 


5.533 


6.794 


8.0567 


9.319 


M 2 tr (asym) a 








6.754 


7.891 


9.026 



a ) see the footnote to Table [Tl 



M{nL) = M (nL) + A str (nL) (L < 4), (6.5) 

then the squared mass 



4(7 3 / 2 L(L + l)(l/p) 
M\L) = M 2 (L) 1 Mq A ' P1 + A 2 t , (6.6) 

If one neglects A 2 tr in Eq. ( |6.6|) which is small (< 0.016 GeV 2 for L < 4) and uses the 
approximation (|6.1| ) for M5 2 (L) then for the orbital excitations with n = the squared mass 
Eq. ( |6.6|) becomes 

M 2 (L) = 8aL - a ^ yP^ L + 1 l + 37r(T = (c^)- 1 L + 37ra (6.7) 

'L + 3tt/8 



where the inverse Regge slope in Eq. ( |6.7| ) is 



(a' L ) 1 = I 8 - v "V /r/v = ^ I a = (6.95 ± 0.02)a. (6.8) 
>L + 3n/8 



y/2(l/ P )(L + l) 



The values of (ct^) 1 are practically constant, see the numbers in Table [TX| , varying from 
the value 6.930 tr for L = 1 to 6.970 tr for L = 4 and we take here (O -1 = (6.95 ± 0.02) a. 
Then 

M 2 (L) = 6.95aL + 3na (6.9) 



or 

f2 



L = 0.UAM 2 (L)/a - 1.358. (6.10) 
It gives for a = 0.18 GeV 2 the Regge slope 



[a 



1.25 GeV 2 or a' L = 0.80 GeV -2 , (6.11) 



in good agreement with the experimental number given in Eq. (|2.3j ) a^(exp) = 0.81 ± 0.01 
GeV -2 . Thus, due to the string corrections we have obtained the correct Regge slope for 
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the spin-averaged masses. However, the intercept in Eq. ( |6.10|) has a very large magnitude 
and an additional contribution to the meson mass must be taken into account. We discuss 
this contribution in Sect. |V11| 



B. Case B. Large L 

For large L the extremal value of the operator v is not equal to or but turns out to 
depend on the parameter (3 as well as on the operator /J,(t). In this case it is a difficult 



problem to find the exact eigenvalues M(asym) of the Hamiltonian Hr, therefore in Ref. |TT 
the eigenvalues of have been calculated in the quasiclassical approximation with the 
following result, 



M 2 (asym) = 2naJL(L + 1) + 3ttc7. (6.12) 



Here, in the asymptotic mass formula (|6.12|) the string correction is already taken into 



account and the constant 3ira is kept to match the solutions for large L to those for L < 4. 
Now, for comparison one can formally define the string correction for large L as the difference 
between the asymptotic mass Eq. ( |6.12|) and the unperturbed mass Mo(nL) Eq. ( |5.2|) 



A str (asym) = \ 3na + 2tiW-L(L + 1) - M Q (L), (L»l,n = 0). (6.13) 



The asymptotic masses are less than M (L) Eq. (|4.6|) for L > 4. The magnitude of A str 
is increasing with growing L and for L = 6 A str (asym) is already ~ 220 MeV. 

From the numerical values of A str (asym) (see Table |m| ) one can see that for L = 4 both 
string corrections, from the asymptotic formula Eq. ( |6.13D and from Eq. ( |6.5[ ), practically 
coincide and in what follows the string correction will be taken in the form ( |6.5|) for L < 4 
and from Eq. fl6.13|) for L > 5 ( when the masses M(asym) are smaller, see Table p| ). 



For L » 1 the Regge slope in Eq. |U3|) is ^Tra)" 1 , i.e. for a = 0.18 GeV 2 , a' L (L » 
1) = 0.88 GeV -2 is larger than for L < 4 and coincides with a'j for the p-tracjectory. Such 
a picture is partly seen in experiment, where for L = 5 the difference M 2 (a§) — M 2 (ps) is 
relatively small and corresponds to the large value a' L ~ 1.1 GeV -2 . However, this growth 
of a' L is likely to be connected with another reason-an effective decreasing of the string 
tension at large distances due to new channels being opened. This effect is considered in 
our paper B. 
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The calculated meson masses (see Table [TTT| ) still are large compared to experiment and 



to get agreement between them a negative constant (a fitting parameter) must be added to 
the squared mass M 2 (nL) |13[. Here we shall not introduce a fitting constant, but instead 
take into account the quark self-energy correction to the meson mass. 



VII. THE QUARK SELF-ENERGY CONTRIBUTION AND MESON MASSES 

Recently it was observed that a negative constant must added to the meson mass, which 
comes from the nonperturbative quark self-energy contribution created by the color magnetic 



moment of the quark [|14}| . This constant is rather large and was calculated with the use 
of the Fock-Feynman-Schwinger representation of the quark Green's function. The total 
nonperturbative self-energy contribution, both from the quark and the antiquark, was found 
to be fully determined by the string tension and by the current mass (flavor) of the quark: 

A SE (nL) = -— -. 7.1 

TT/j, (nL) 

Here no(nL) is just the constituent mass defined by Eq. fl5.1|). The constant r/(/) depends 
on the flavor: its numerical value for a quark of arbitrary flavor was calculated in Ref. [|14| , 
in particular for the light mesons we take as in Ref. [|H| 



v (nn) = 0.90. (7.2) 



The self-energy terms, as well as the meson masses, are given in Table [TV] for the ground 
states (n = 0, L < 5) from which one can see that Ase(^) decreases as a function of n and 
L, being proportional to fiQ l {nL). Still it is rather large (equal to -300 MeV) even for L = 5. 

With the self-energy and the string corrections taken into account the spin- averaged 
meson mass M(nL) is fully determined. The Coulomb correction will be discussed in the 
next Section and calculated in Appendix [FJ. 

The meson mass is now given by 

M { nL) = M„(nL) - ~^,(L< 4) (7.3) 

and for the linear potential can be written as 

»->/ »,/ 2a(l/r)L(L + l) War) frr A . 

M(nL) = A M nL) - > - - J, (7.4) 
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TABLE IV: The nonperturbative quark self-energy correction Ase(-^) and the meson masses M(L) 
in GeV for the ground states (a = 0.18 GeV 2 , r\ = 0.9). 



L 





1 


2 


3 


4 


5 


Ase(L) 


-0.616 


-0.460 


-0.383 


-0.335 


-0.301 


-0.294 


M(0,L) 


0.723 


1.279 


1.685 


2.017 


2.30 


2.514 a 


Mexp(0,L) 


0.612 


laj(1.252) 


7T 2 (1.67) 


04(2.014) 


P5 


a 6 (2.45±0.13) 






l/j(l-245) 


103(1.69) 
ws(1.67) 


/ 4 (2.034) 


2.33 ±0.04 


/ 6 (2.47± 0.050) 



) this mass was calculated from the asymptotic formula (|6.1 



using the relations (|5.2|) and ( |5.8| ). The calculated meson masses (L < 4) coincide with good 
accuracy with the experimental values ( see Table |V| ). 
For large L (n = 0) 

M — J 37rcr + 2irayjL(L + 1) + A SE (L). (7.5) 



VIII. THE INTERCEPT OF THE REGGE TRAJECTORY 



From the mass formula ( [7.4|) it follows that the self-energy term enters M(nL) in such a 
way that the negative constant Co, 



C 



TV 



appears in the squared spin-averaged mass M 2 (L): 



M 2 (nL) = (M + A s 



^ 2 32ar/ t /16cttA 2 
V^mJ ' 



(8.1) 



(8.2) 



7T \ /1 1V1Q 

Here the terms A 2 tr and A str Ase will be neglected, because they give small contributions 
for L < 4, while the term A| E is kept, since it is not small in all states. The constant Co 
is rather large and for a = 0.18 GeV 2 is equal to —1.65 GeV. Using the expression ( |6.6| ) for 
the mass (M + A str ) 2 and Eq. ( |B~7Ll) for M5 2 , Eq. ( ^.2[ ) can be presented as 



M 2 (L) = ( a ' L y l L + b(L), 



(8.3) 



with 



b(L) = a 



3vr 



32r/ 2 



32r/ 

~V + 7T 2 (L + 37T/8) 
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b(L = 0) = a 



32r? 256r? 2 

37T h 



(8.4) 



7T 7T 2 M 2 (1S)_ 

where for M (L = 0) it is better to use the exact value, M (1S) = 3.157 -y/cr and a' L was 



already defined by the expression (|6.8|) . From (|8.3|) the intercept is 



a L (0) = a L (M 2 = 0) = -a' L b(L = 0) = ° } . (8.5) 

o.95o~ 

Note that in b(L) the combination (37r — 32i]/ir)a is a small number (equal 0.046 GeV 2 for 
a = 0.18 GeV 2 ) and therefore for the intercept the contribution of the self-energy term A| E 
is dominant. 

From Eq. ( |3.4| ) it is clear that b(L) is sensitive to the value of the flavor factor rj, which 
may introduce an uncertainty on the order of 5%. 

With the use of the analytical expression flS.4| ) and the exact value of M (L = 0), t](nn) = 
0.90 the quantity b{L = 0) is equal to 

b(L = 0) = 2.365(7. (8.6) 

Then the intercept given by Eq. ( p75|) takes the value 

a L (0) = -(a' L ) b(L = 0) = -2.356/6.95 = -0.34. (8.7) 

This number is in good agreement-larger by 10% only- with the experimental value 
° ; l(0) = —0.30 ± 0.02. It is essential that the intercept does not depend on the string 
tension but instead is sensitive to the flavor parameter rj. Just for this reason the intercept 
for the mesons with different flavor depends on the flavor. 

So, finally, the Regge L-trajectory calculated in the QCD string approach with a = 0.18 
GeV 2 is fully determined, 

L = 0.80 M 2 (L) -0.34 (8.8) 

and appears to be very close to Eq. ( ^3|) obtained from a fit to the experimental spin- 
averaged meson masses, see Fig. H From Eq. the averaged mass M(n — p) is found: 

M 2 (1S) = 0.425 GeV 2 or M (IS) = 0.652 GeV, (8.9) 

which corresponds to a 7r-meson mass M(ir) = 301 MeV. This number turns out to be 
smaller than M(1S) = 0.723 GeV calculated directly from Eq. ( |6.5|) and this discrepancy 
illustrates how sensitive M(1S) is to the approximations used. 
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IX. COULOMB INTERACTION 



In the previous sections good agreement of the spin-averaged meson masses (for the 
ground states with L / 0) was obtained without taking into account the Coulomb interac- 
tion. It is of interest to check whether the Coulomb effects are actually suppressed for L > 
states and how large is Coulomb correction to M(tt — p). 

To this end we solve the Salpeter equation with the string potential taken as a linear plus 
Coulomb term, i.e., with the Cornell potential: 

Va{r) = ar- (9.1) 
3 r 

where a s = constant can be used, since the light mesons have very large sizes, R > 1.0 fm, 



and at such distances the strong coupling is saturated and close to the "freezing" value |L5 
If for the string tension one takes a = 0.19 GeV 2 , then the fitted value of a s appears to be 



just the same as for heavy quarkoniapO|, |21 



a s = 0.39 (9.2) 

However, the masses of the ground states, including the IS" state, can be nicely described 
with a smaller value for the coupling constant, 0.20 < a s < 0.39, if correspondingly the 
value of a is taken from the range 0.18 GeV 2 < o < 0.19 GeV 2 . 

The main characteristics of the qq system like the eigenvalues Mc{nL) of Eq.( [4.6| ) using 
the Cornell potential, the constituent masses fic( n L) defined by Eq. ( |5.1| ) together with the 
string and the self-energy corrections are presented in Appendix [B] in Tables [X| and BO. Here 
in Table we give only the results of our calculations for the spin-averaged masses Mc(nL). 
Note that in the Coulomb case the relation (|5.2j ) is not valid and therefore the meson mass 
Mc{nL) as well as A str and Ase should be written through the constituent mass (denoted 
as ncijiL)) as in Eq. Q7.3Q (see Table |X| where the eigenvalues are given for a = 0.19 GeV 2 , 
77 = 0.90, and a s = 0.39). 



With the use of the string and the self-energy corrections from Table [XI] the spin-averaged 
meson masses Mc(L), Eq. (|7.3j ), are determined and their values are given in Table 
together with the experimental numbers. 

If now one compares the meson masses Mq{L) with those for the linear potential from 
Table PH[ then one can see that in the Coulomb case for the 15* and IP states a better 
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TABLE V: The spin-averaged masses Mq{L) in GeV, theoretical and experimental, for the ground 
states (n=0) (a = 0.19 GeV 2 , a s = 0.39). 



Jj 


n 


1 


2 


(J 


4 




M C {L) 


0.632 


1.220 


1.650 


2.00 


2.29 


2.51 


M exp (L) 


0.612 


M(fj) = 1.24 


vr 2 (1.66) 


04(2.014) 


p 3 (2.30) 


a 6 (2.45) 








p 3 (l-69) 


/ 4 (2.03) 




/e(2.47) 






M(aj) = 1.25 


w(1.65) 














w 3 (l-67) 









TABLE VI: The exact and perturbative Coulomb corrections Ec(L) (in MeV) 



L 





1 


2 


3 


4 


5 


£"(7 (exact) 


-219 


-132 


- 103 


- 86.7 


- 76.2 


- 68.8 


^c(pert) 


-194 


-126 


-99.4 


- 84.9 


- 75.1 


- 68.1 


difference 


11.4% 


4.5% 


3.5% 


2.0% 


1.4% 


1.0% 



agreement with the experimental numbers is obtained, however, in the Coulomb case the 
string tension appears to be larger, a = 0.19 GeV 2 . The calculated mass Mc{lS) = 0.632 
GeV is very close to the value Eq. (|8.9|) from the Regge trajectory Eq. ( |8 . 7| ) . 

Now the Coulomb correction can be formally defined as the difference between the exact 
eigenvalues, Mc(L) and M(L): 



^c(exact) = M c (nL) - M(nL) 



and compared with the Coulomb corrections E c (pert): 



£ c (pert) = -~a s (l/r) 



(9.3) 



(9.4) 



obtained when the Coulomb interaction is considered as a perturbation (see Table jVT|) . In 
Eq. ( |9.4|) the m.e. < 1/r > is to be taken for the linear potential with the same a as in the 
Cornell potential. 



The numbers in Table demonstrate that the exact and perturbative corrections coin- 
cide with an accuracy better than 5% for all states with L > (for the IS" state the difference 
is 11% ) and therefore these corrections can be calculated as a perturbation. 
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For the nL states one should also take into account the difference between the exact 
constituent mass Hc{L) and fio(L) for the linear potential; they are related as follows 



Hc{lL) « fi (lL) + \E c \/3 (n = 0), 

fic(nL) « Vo{nL) + \E c \/4, (n^O). (9.5) 
This correction to the constituent mass is mostly important for the IS state. For larger 



n the difference between fic and /xo can be neglected. As seen from Table [VT], due to the 
Coulomb interaction all masses are shifted down by an amount in the range of 70 to 200 
MeV and therefore a larger value of a is needed, a = 0.19 GeV 2 for a s = 0.39, than for the 
linear potential. 

However, one cannot take an arbitrary or too large value for a, otherwise the Regge 
slope a' L would be small and in contradiction with the experimental value. Therefore, in the 
Coulomb case only values a = 0.19 ±0.10 GeV 2 are allowed. Then to obtain agreement with 
experiment using a < 0.20 GeV 2 a restriction on the value of the strong coupling constant 
is found: 

a s < 0.40 (a < 0.20 GeV 2 ), (9.6) 

otherwise correct numbers for the Regge slope and the intercept cannot be obtained simul- 
taneously. 

This upper limit ( |9.6|) for a s appears to be in accord with the freezing value of the two- 
loop a B {q 2 = 0) = 0.45 (with the QCD constant A® = 330 MeV, N f = 3) obtained in 
background field theory fT5|j . 



X. CONCLUSIONS 



In the framework of the QCD string approach the spin-averaged meson masses with 
L < 5 (n = 0) have been calculated and expressed through a single parameter-the string 
tension a and a set of universal numbers. In this approach the kinetic energy is of the same 
type as in the spinless Salpeter equation. The constituent mass and the nonperturbative 
quark self-energy are calculable and also depend on the string tension only. 

This is the first time accurate predictions for the meson masses have been obtained relying 
on one parameter only, that is directly connected to the confinment mechanism in QCD. 
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The analytical expressions for the slope and the intercept of the Regge L-trajectory 
(when the spin splittings are not taken into account) have been deduced, giving rise to a 
value a' L = (6.95cr) _1 = 0.80 GeV -2 (a = 0.18 GeV 2 ) which coincides with the experimental 
number. This L-trajectory can be considered as a universal one since in the approximation 
of closed channels it does not depend on spin and isospin . 

It is shown that the Regge intercept does not depend on a and a(M = 0) (theory) = —0.34 
turned out to be only 10% larger than a(M = 0)(exp) = —0.30 ± 0.02. From this intercept 
M(1S) = 652 MeV corresponds to a 7r-meson mass equal to 300 MeV (chiral effects have 
been neglected here). 

For all orbital excitations with L ^ the calculated masses are in a good agreement with 
existing experimental data. 

In order to obtain this good agreement with the data we find it necessary to impose a 
restriction on the value of a s that is in accord with the freezing picture. 



The authors would like to express their gratitude to the theory group of Thomas Jefferson 
National Acelerator Facility (TJNAF) for their hospitality. This work was partly supported 
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APPENDIX A: DETAILED SPECTRA 



The eigenvalues and the wave functions of the SSE equation were calculated with the help 
of the code used before j^T], [22 . The eigenvalues and relevant matrix elements are given in 



Tables [VI I| and [VIII| for the linear potential and in Tables PQ-|XI| for the Cornell potential. 

From Table [Villi one can see that the difference between the m.e. < 1/r 3 > for the SSE 
and the Airy equations for the P-wave states turn out to be large reaching 40% for n > 2. 



As briefly discussed in Sect. Ill the reason behind these differences lies in the different 



asymptotic behaviors of the eigenfunctions of these two equations. In Table |IX| we also give 
the constituent masses Eq. ( |4.6| ) and the m.e. < 1/r > and < 1/r 3 > for the ground states 
(n = 0) with L < 6. 

The calculated m.e. < 1/r > is used to obtain the string and Coulomb corrections, while 
the m.e. < 1/r 3 > can be used to calculate the hyperfine and fine-sructure splittings for the 
mesons with L ^ 0. 
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TABLE VII: The eigenvalues Mo(nS), constituent masses no(nS), p,a(nS), and matrix elements 
< 1/r > (in GeV ) for the Salpeter Eq. ([4.6|)and the Airy equation Eq. ( |5.4| ) with the linear 
potential ar (a = 0.18 GeV 2 , L = 0). 



n 





1 


2 


3 


4 


M (nS) 


1.3394 


1.9980 


2.4985 


2.9151 


3.2797 


Ho(nS) 


0.3348 


0.4995 


0.6246 


0.7289 


0.8199 


fio(nS) 


0.3519 


0.5351 


0.6703 


0.7826 


0.8807(3) 


< 1/r > (SSE) 


0.3638 


0.3299 


0.2959 


0.2734 


0.2559(5) 


< 1/r > a (EA) 


0.3328 


0.2669 


0.2334 


0.2118 


0.1996(5) 



a ) These m.e. are calculated from the Eq.(18) with j~Lo(nS) definec 



by Eq.(20). 



TABLE VIII: The matrix elements < 1/r 3 > (GeV 3 ), mass eigenvalues M (nP) (GeV), and 
constituent masses hq{uP) and p,o(nP) (GeV) for the P-wave states (a = 0.18 GeV 2 ). 



n 





1 


2 


3 


4 


M (nP) 




1.7924 


2.3153 


2.7505 


3.1291 


3.4682 


MoK) 




0.4481 


0.5788 


0.6876 


0.7823 


0.8671 


fio(nP) 




0.4620 


0.6115 


0.7320 


0.8335 


0.9278 


< 1/r 3 > (SSE) 




0.0264 


0.0422 


0.0539 


0.0635 


0.0718 


< 1/r 3 > (EA)= 


/iff/4 


0.0208 


0.0275 


0.0329 


0.0376 


0.0417 



a ) see the footnote to Table |VT 



From Table [X| one can see that in the Coulomb case the constituent masses //c(lS') and 
/ic(l-P) are larger by 29% and 10% respectively, than for the linear potential (see Table |V111[ ) 
and therefore Ase is smaller for them. 
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TABLE IX: The constituent masses (GeV) and the matrix elements < l/r k > (GeV fe ), (k=l,3), 
of the SSE for the ground states(n=0) (a = 0.18 GeV 2 , L < 6). 



L 


1 


2 


3 


4 


5 


6 


W)(0L) 


0.4481 


0.5387 


0.6163 


0.6853 


0.7480 


0.8058 


< 1/r > 


0.2362 


0.1867 


0.1589 


0.1406 


0.1274 


0.1173 


< l/p > Vl + 1 


0.787 


0.762 


0.742 


0.741 


0.736 


0.732 


< 1/r 3 > 


0.0264 


0.0098 


0.0054 


0.0035 


0.0026 


0.0019 



APPENDIX B: RESULTS FOR THE CORNELL POTENTIAL 

In this appendix we present some auxiliary values for the Cornell potential. 

TABLE X: The eigenvalues Mc(L), the constituent masses /Uc*(L) and the m.e. < 1/r > (in GeV) 
for the SSE with the Cornell potential for the ground states (n = 0). (cr = 0.19 GeV 2 , a s = 0.39) 



L 





1 


2 


3 


4 


5 


M C (L) 


1.157 


1.710 


2.111 


2 446 


2.740 


3.005 


»c(L) 


0.415 


0.496 


0.580 


0.656 


0.710 


0.745 


< 1/r > L 


0.484 


0.266 


0.202 


0.170 


0.145 


0.130 



TABLE XI: The string and self-energy corrections (in GeV) for the SSE with the Cornell potential 
(a = 0.19 GeV 2 , a s = 0.39). 



L 





1 


2 


3 


4 


5 


Astr 





-0.051 


-0.086 


-0.112 


-0.075 


-0.068 


A SE 


-0.525 


-0.439 


-0.375 


-0.332 


-0.307 


-0.282 
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